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We show numerically that, in a Josephson ladder with periodic boundary conditions and subject 
to a suitable transverse magnetic field, a vortex excitation can spontaneously break up into two 
or more fractional excitations. If the ladder has N plaquettes, and N is divisible by an integer q, 
then in an applied transverse field of 1/q flux quanta per plaquette the ground state is a regular 
pattern of one fluxon every q plaquettes. When one additional fluxon is added to the ladder, it 
breaks up into q fractional fluxons, each carrying 1/q units of vorticity. The fractional fluxons are 
basically walls between different domains of the ground state of the underlying 1/q lattice. The 
fractional fluxons are all depinned at the same applied current and move as a unit. For certain 
applied fields and ladder lengths, we show that there are isolated fractional fluxons. It is shown 
that the fractional fluxons would produce a time-averaged voltage related in a characteristic way to 
the ac voltage frequency. 

PACS numbers: 74.81. Fa,74.50.+r,74.25.Qt 



I. INTRODUCTION 

Ladder arrays of Josephson junctions have been exten- 
sively studied both theoretically and experimentally (see, 
for example, Refs. 0, H H H Q, @, H El El El 111 El 
E[HiE3) They are of interest, in part, because they 
show a rich variety of equilibrium and dynamical behav- 
ior. In addition, they are an important testing ground 
for concepts of phase transitions and quantum behavior 
in lower dimensionality. The interest in such materials is 
even greater now, because advances in microfabrication 
allow such ladders to be made almost to specifications, in 
order to test theoretical predictions or possibly to make 
new types of Josephson devices. 

Recently Chandran and Kulkarni[lj| considered the 
behavior of a hypothetical Josephson ladder array with 
alternating 7r and junctions. A n junction is one in 
which the Josephson supercurrent / = 7 c sin(A0 + it), 
where A(f> is the phase difference across the junction 
and I c > 0. Such junctions can be formed in a vari- 
ety of ways, e.g., by controlling the energy distribution 
of the current-carrying states in the normal metal within 
the iunction|l9|. by preparing a junction between two 
cuprate superconductors in the presence of bound states 
in the interface materially, or by connecting two su- 
perconductors across a ferromagnetic layer [2l| . Thus, an 
array of alternating and 7r junctions could, in principle, 
be made in the laboratory. In Ref . E3 > it was shown that 
when a 2tt fluxon is introduced into such a — it ladder, 
it will break up into two separate tt fluxons, each of mag- 
nitude $o/2- This fractionalized vortex was predicted to 
have unusual current-voltage characteristics which could 
readily be detected experimentally. 

In this paper, we show that similar fractionalized vor- 
tices are produced in conventional ladder Josephson ar- 
rays made of junctions, in a transverse magnetic field. 
Specifically, we consider a ladder having N square pla- 
quettes, where N is divisible by an integer q, in a trans- 
verse magnetic field. If that field is of magnitude equal 



to one flux quantum per q plaquettes, the ground state 
of the ladder is a periodic array of N/q equally spaced 
fluxons. If now one fluxon is added to that ladder, we 
find that the added fluxon breaks up into q fractional 
fluxons, each of magnitude &o/q. If q = 2, the fluxon 
pattern is similar to that found in Ref. E3j but we find 
analogous patterns for all other values of q which we have 
tested. These extra fluxons have IV characteristics and 
vorticity patterns which should be readily detectable ex- 
perimentally. Furthermore, the fractionalized fluxons are 
typically "confined" that is, they are typically depinned 
at the same current and move at the same velocity. For 
an array in which N is not divisible by an integer q, we 
show numerically that it is possible to produce isolated 
fractional fluxons, by applying a transverse magnetic field 
of suitable magnitude. 

The remainder of this paper is arranged as follows. In 
Section II, we develop the equations of motion from a 
standard Lagrangian for the Josephson ladder. In Sec- 
tion III, we give our numerical results, which show that 
fractional vortices are generated in the ladder at suit- 
able applied transverse fields. Finally, in Section IV, we 
briefly compare our results to other models which give 
rise to fractional excitations, and give a concluding dis- 
cussion. 



II. FORMALISM 

We consider a Josephson ladder with periodic bound- 
ary conditions, as shown schematically in Fig. The 
ladder is assumed to consist of a collection of small junc- 
tions, each of critical current I c , which are inductively 
coupled together. The edges of the ladder are parallel to 
the x axis, while the rungs are in the y direction. A dc 
current / is injected into each junction on one edge of 
the ladder and extracted from each junction on the other 
edge. A magnetic field B = Bi is applied perpendicular 
to the ladder. The geometry is readily achievable experi- 
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mentally, in the form of a coplanar ring perpendicular to 
B. 

In the absence of dissipation, it is convenient to de- 
scribe the ladder by the following Lagrangian: 
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Here K represents the charging energy of the small junc- 
tions, U = 4e 2 /C, and C is the capacitance of the small 
junction. J\ = HI c /2e is the Josephson coupling energy 
of the small junction, J 3 = HI/(2e), and J 2 = H 2 /(8e 2 L), 
where L is the inductance of the wire segments coupling 
points (i,£) to (i + 1,£) and (i,u) to (i + l,u). Finally, 
if we use the gauge A = Bxy, where A is the vector 
potential, then A, = (2tt/$ ) J-f A • dl = (2n/<i> Q )a 2 iB, 
where $o = hc/(2e) is the flux quantum, a is the lattice 
constant (see Fig. 1), and the subscripts I and u indicate 
the lower and upper ends of the junction, as shown in the 
Figure. The dots indicate time derivatives. Note that we 
neglect the self-induced magnetic field, by taking B equal 
to the externally applied field. 

If we now make the substitutions fa = 9 i u — 9a, Xi — 
(@i,e + 0i,u)@>E|) we ma y reexpress the Lagrangian as 
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The Lagrange equations of motion are then 
(d/dt(d£/d(j)i)) — d£/d(f), with an analogous equa- 
tion for the variables Xi- Since £ is independent of Xii 
the equation of motion for Xi leads to the condition that 
Xi — Xi+i is independent of i. After some algebra, it is 
readily shown that the equation of motion for (pi is 

(j>i = - sin(0j + 2nif) + \j(<f>i-i - 2<fii + cf) l+ \) 

+ {l/i c )4>i-j>i/Qj- (5) 

Here, we have introduced the frustration / = £?a 2 /$o, 
and a dimensionless time r = u) p t, where u p = 
[2e/ c /(/iC)] 1 / 2 is the junction plasma frequency, and 
Xj = (J21 1 J1) 1 / 2 is a dimensionless Josephson penetra- 
tion depth. The dot now represents a derivative with 
respect to r. Finally, we have added an additional 
dissipative term —4>i/Qj by hand to the equation of 



FIG. 1: Schematic of the ladder modeled in this paper. We 
show the i — 1, i, and i + 1 rungs, u and I denote the upper 
and lower edges of the ladder. An external magnetic field is 
applied in the z direction; a dc driving current / is applied to 
the lower end of each rung in the y direction and extracted 
from the upper end. The ladder is assumed to have the topol- 
ogy of a ring (periodic boundary conditions). The full lines 
along the upper and lower edges of the ladder (parallel to x) 
are treated simply as inductors in our model. 



motion, where Qj is a dimensionless junction quality 
factor, assumed to be the same for all junctions. In 
the resistively and capacitively shunted junction model 
(RCS.Ji 22 . Qj = (2eR 2 I c C/h) 1 / 2 , where R is the junc- 
tion shunt resistance. 



III. NUMERICAL RESULTS 

A. Numerical Method 

In most of our calculations, we have studied a ladder 
with N = 120 rungs and periodic boundary conditions. 
This number of rungs was chosen because it is the small- 
est integer divisible by 2, 3, 4, 5 and 6, and is sufficiently 
large to allow the fractionalization to be apparent in the 
simulations. We solved our set of 240 first order, non- 
linear, differential equations using a constant-time-step, 
fourth-order Runge-Kutta method |2^|. We have gener- 
ally considered underdamped junctions, with Qj = 10, 
and values of A 2 between 0.01 and 10. We initialized the 
system in a random phase configuration at r = and 
I / I c = 0. The equations are integrated forward in time 
up to r eq = 3000, by which time we assume the system 
has reached a steady state. We then calculate the time- 
averaged voltage drop across a rung of the ladder (also 
averaged over the rungs) by averaging over an additional 
dimensionless time units, i. e., 
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= 5000. The dc driving current is then in- 
creased by 0.01J C and the calculation is repeated. We 
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FIG. 2: (a) IV curves for / = 1/2, Qj = 10.0, and Aj = 
1.0 for increasing and decreasing current (the directions are 
indicated by arrows). There are two voltage jumps in the 
IV curve on the increasing branch: at the first, the fluxon 
lattice is depinned; at the second, the junctions switch to the 
voltage state, (b) Schematic of ground state vortex number 
configuration for the f — 1/2 ladder at I/I c = 0; the +'s and 
blank squares represent vortex numbers n = 1 and n = 0. 
Only 40 plaquettes are shown. 



continue ramping up the current by steps of 0.01 until 
I/I c = 1.2, then decrease I/I c in steps of 0.01 back to 
zero. These calculations give the IV curves shown below 
in Figs. Efa), El Hand El 



1/2 



In Fig. Efa), we show the full IV curve for / = 1/2, 
Qj = 10, and X 2 j = 1. The directions of the arrows 
indicate whether the current is increasing or decreasing. 
The IV curve is clearly hysteretic for this value of Qj, and 
has two discontinuous jumps on the increasing current 
branch. At the lower jump, near I/I c — 0.13, the system 
jumps into a state where the fluxon lattice is depinned 
and starts to move through the ladder as a unit, giving 
rise to a finite voltage across the ladder. At the upper 
jump, near I / I c = 0.5, all the junctions switch to a finite 
voltage state and the fluxon excitations are expelled from 
the ladder. 

To represent the fluxon lattice pictorially, we use the 
concept of a vortex number\2^. The vortex number of 



the a plaquette is defined as 
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Here the gauge-invariant phase difference for each leg of 
the plaquette is written (9i — 9j —Ay) and is defined to lie 



in the range [— tt, it]; the summation is taken around the 
plaquette in the counterclockwise direction (viewed from 
the positive z axis). Thus n a is or a positive or negative 
integer for each plaquette. In the lowest energy state, we 
expect that ^ Q n a — Nf, where N is the number of 
plaquettes, and also that all the individual n a = or 1. 

The two values of n a can generally be distinguished by 
the current pattern around the plaquette: n a = 1 and 
correspond to a clockwise and counterclockwise current 
flow, respectively. The expected ground state for / = 1/2 
at I/I c = is shown pictorially in Fig. Etb). The +'s 
and empty squares represent plaquettes of n a = 1 and 0, 
respectively. 

In practice, we have not found a way to calculate n a 
directly, because we do not compute the variables Xi 
mentioned above. Therefore, in order to confirm the 
vortex pattern at various fields and currents, we cal- 
culate the "normalized flux," defined by the relation 
27r < I>i/$o = 4>i ~ 4>i+i + 2tt /, in each plaquette. We show 
a number of examples of this normalized flux below. 

When the vortex lattice is depinned near I/I c = 0.13, 
the entire lattice of +'s moves as a unit through the lad- 
der, generating a voltage. In Fig. El we show the space- 
averaged voltage 
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across the ladder for 4500 < r < 4520, for the parameters 
of Fig. Eta) at I/I c = 0.13, just above the lower jump in 
the IV characteristics. We interpret the periodic voltage 
oscillations as arising from the motion of the fluxon lat- 
tice through the bumpy "egg-carton" potential provided 
by the plaquettes of the Josephson ladder: the velocity 
of the fluxon lattice varies periodically in time because 
the fluxon lattice moves alternately faster and slower as 
it moves through the steeper and less steep part of the 
egg carton p5|. 

To confirm this picture, we note from Fig. El that the 
dimensionless period of the voltage oscillation, which we 
denote At, is ~ 2.7. We expect that this should be the 
dimensionless time required for the fluxon lattice to ad- 
vance one lattice constant a. Now if the fluxon lattice 
moves a distance a during this period, then the aver- 
age speed of the fluxon lattice is v = (Na)/T, where T 
is the time required for the lattice to complete one cir- 
cuit around the ladder. During this time, each junction 
should experience a phase change of 2ir N /, since Nf 
fluxons cross each junction during one circuit. Thus, the 
time-averaged voltage (V) across each junction should 
be [h/(2e)]2%Nf /T. Introducing a dimensionless voltage 
(V) / (RI C ), and noting that the time required to move one 
lattice constant a is uj p T /N — At, we obtain 
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Using the parameters of Fig. 
I/J c = 0.13 and At = 2.7, eq. 



Eta), together with 
© gives (V)/(RI C ) = 
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FIG. 3: Time-dependent voltage V(t)/(RI c ) across one rung 
of the ladder, averaged over the rungs, and plotted for dimen- 
sionless time r between 4500 and 4520, for / = 1/2 and the 
parameters of Fig. HJa). The dimensionless period of voltage 
oscillation is ~ 2.7 at I/I c = 0.13. 
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FIG. 4: Normalized "flux" 27r$;/$ (as defined in the text) 
vs. plaquette number i, for the parameters of Fig. [^Ja) and 
I/I c = 0. For clarity, only 60 of the 120 plaquettes are shown. 



2tt(1/2)/(10 • 2.7) = 0.116. The actual time-averaged 
voltage read off from Fig. Efa) at I/I c = 0.13 is 0.117, in 
excellent agreement with this calculated value. Thus, the 
time-dependent voltage shown in Fig. ©is consistent with 
the picture of a fluxon lattice moving as a unit through 
the washboard potential provided by the ladder. For cur- 
rents slightly higher than I/I c = 0.13, this picture of a 
moving vortex lattice continues to hold. As the current 
is increased still further, for Qj — 10, the voltage is no 
longer perfectly periodic in time (not shown in the Fig- 
ure). For I/I c greater than about 0.51, the individual 
junctions switch to a voltage state and the moving fluxon 
lattice disappears. 

In Fig. © we plot the calculated "normalized flux" 
$ 4 /$ for f = 1/2 and I/I c = 0. Clearly, $,/$ alter- 
nates between two values, indicating that the flux lattice 



is periodic with a period of two unit cells, as suggested 
by Fig. 2(b). The two ground state fluxes are approxi- 
mately 27r/±/ with / = 1/2. A similar pattern has been 
seen previously in Ref. Jig for a ladder with alternating 
and 7r junctions, and zero applied magnetic field. 



C. f = 1/2 +1/N 

Next, we consider the case / = 1/2 + 1/120, equivalent 
to adding one fluxon to the f — 1/2 ladder. The calcu- 
lated results are shown in Figs. [5] and© Fig. ©shows 
the IV characteristic as computed on increasing and de- 
creasing the current in steps of 0.01/ c . On the increasing 
branch, there are now three steps in the IV curve. Two of 
the steps are easily seen from the Figure, but the lowest 
jump, which occurs near I/I c = 0.01, is less visible and 
shown in the inset for < I/I c < 0.15. The IV curve is 
hysteretic in two regions of current, but non-hysteretic at 
the lowest non-zero voltages. On the decreasing branch 
there are still three distinct regimes of non-zero voltage. 
We interpret these as follows: (i) at the lowest currents, 
a single fluxon is depinned and moves through the fluxon 
lattice formed by the remaining 60 fluxons; (ii) at higher 
currents, the fluxon lattice is depinned, as at / = 1/2, 
generating a larger voltage; and finally (iii) at still higher 
voltage, all the junctions switches to a voltage state and 
the ladder jumps up to the resistive branch. 

Fig. ©a) shows the flux pattern 27r<I> J ;/$o for / = 
61/120 and I/I c = 0, keeping the rest of the parameters 
the same as in Fig. Eta)- The f — 1/2 pattern is now 
distorted in two regimes, near i = 50 and i = 110, while in 
the rest of the ladder the pattern is similar to that of / = 
1/2. To compare the two patterns more quantitatively, 
we note that the sums $i/$o for Figs. ©and©a) are 
60(27r) and 61(27r) respectively. This result shows that 
the effect of the additional magnetic field is to add one 
additional fluxon to the ladder. However, this fluxon does 
not enter the ladder as a unit. Instead, it is fractionalized 
into two 1/2-fluxons, each of which carries 1/2 a unit of 
flux. These two 1/2-fluxons obviously repel one another, 
since they prefer to be situated as far away from one 
another as possible, i. e. about 60 plaquettes apart in 
the 120-plaquette ladder. 

To suggest intuitively why the extra fluxon fraction- 
alizes in the ladder, we show in Figs. © (b) and (c) two 
possible scenarios for how the extra fluxon enters the lad- 
der. In (b), the fluxon enters as a unit, without distorting 
the underlying f = 1/2 fluxon lattice. But this scenario 
leads to n = 1 vortices in three consecutive plaquettes, 
as indicated by the three adjacent plus signs. This is 
obviously a high energy state, since vortices of like signs 
repel one another. A more plausible state is shown in Fig. 
© (c) , where the extra fluxon has fractionalized into two 
1/2-fluxons. In this case, there are only two points in the 
ladder where there are two adjacent + signs, and none 
with three. Our numerical simulations indicate that this 
is indeed a lower energy state than the un-fractionalized 
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FIG. 5: IV curve for / = 61/120, Qj = 10 and X 2 j = 1 for 
increasing and decreasing current (as indicated by arrows). 
The inset shows the IV curve for < I/I c < 0.15. A finite 
voltage in this regime is produced by the two moving 1/2- 
fluxons. 



fluxon shown in Fig. (b). We have also found that 
both halves of the fractionalized fluxon are depinned at 
the same I /I c , and move together. 

Finally, we consider the expected time-dependent volt- 
age for the two moving 1/2-fluxons in the 61/120 ladder. 
As explained above, each 1/2 fluxon is really a kind of do- 
main wall, consisting of two adjacent + vortex charges. 
The expected mechanism by which a 1/2-fluxon moves 
is that one of the two adjacent + charges moves to the 
nearest empty plaquette. This jump causes the domain 
wall to move by two plaquettes. Thus, if the ladder has 
length N plaquettes, a 1/2- fluxon will circulate around 
the ladder twice in N moves. This causes a phase slip of 
2 • 7r, since there is a 7r phase slip each time the 1 /2-fl.uxon 
circulates once around the ladder. Thus, if T is the time 
required for the 1/2-fluxon to circulate once around the 
ladder, the time-averaged voltage is (V) — (h/2e)x/T. 
On the other hand, V(t) for a single 1/2-fluxon should 
have an ac component with a period of \/v = 2T/N, 
where v is the ac frequency, since 2T/N is the time re- 
quired for a 1/2-fluxon to make a single move of two 
plaquettes. 

Combining these two relations, we get (V) = 
\h/(2eN)]v for a single 1/2-fluxon. This is precisely the 
relation between (V) and v which we find for a ladder 
containing a single 1/2-fluxon, as is discussed later in 
this paper. 

For the present case of a 61/120 ladder, there are, as 
already mentioned, two 1/2-fluxons, as shown in Fig. 6. 
Thus, at a current / in the regime where the two 1/2- 
fluxons are moving, we expect that (V) will equal twice 
this value, or 



(V) = 2[h/(2eN)]v, 



(10) 



where v represents the frequency of the ac voltage gener- 
ated by a the motion of a single 1/2-fluxon at the same 
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FIG. 6: (a) Normalized flux vs. plaquette number for the 
parameters of Fig. |21 except that / = 61/120. (b) Schematic 
of an un-fractionalized, high-energy vortex charge configura- 
tion for / = 61/120. (c) Schematic of the actual ground state 
vortex charge configuration at / = 61/120 and I/I c = 0. In 
both (b) and (c), we show only 40 of the 120 plaquettes. (In 
the actual ladder, the two fractional vortices are separated by 
60 plaquettes.) 



current. We can rewrite this equation in terms of Qj and 
At = Up/v with the result 



(V) _ 2 ■ (2tt) 
RI C NQjAt' 



(11) 



In the present case N — 120 so the right-hand side equals 
2tt/(60QjAt). The calculated (V), in the regime where 
the two 1/2-fluxons are moving, agrees very well with the 
prediction of eq. iJUJ. Thus, the IV characteristics in this 
regime are produced by two 1/2-fluxons moving with the 
same average velocity. 



D. 



1/3 



If an added fluxon fractionalizes when added to a lad- 
der at / = 1/2, what happens at other values of /? To 
answer this question, we have carried out similar calcula- 
tions for 1/3, 1/4, 1/5 and 1/6 in an N = 120 ladder, as 
well as for each of these values of / with one fluxon added 
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FIG. 7: IV curve for / = 1/3, Qj = 10.0, A 2 7 = 1.0 and 
TV = 120. The inset shows the IV curve for < I/I c < 0.05. 

or subtracted. The calculated IV characteristics and nor- 
malized flux for / = 1/3 are shown in Figs. |7|and[Sfa). 
The behavior of the IV curves shown in Fig. \7\ is most 
easily understood on the decreasing current branch. For 
the highest currents, the junctions are all in the voltage 
state, and the fluxon lattice has disappeared. At inter- 
mediate currents, the / = 1/3 fluxon lattice is moving 
through the washboard potential provided by the lad- 
der. At the lowest currents, the fluxon lattice is pinned 
and the time-averaged voltage vanishes. On the increas- 
ing current branch, the fluxon lattice appears to depin 
in two stages, corresponding to the two voltage jumps 
visible in the Figure. 

Fig. IBfa) shows the normalized flux $i / $o versus pla- 
quette number i at I/I c = and figure |SJb) shows the 
vortex charge pattern configuration in the equilibrium 
state. The flux pattern is again periodic, as at / = 1/2, 
but now with a period of three plaquettes, with fluxes 
2tt/ + /, 2nf + f, and 2nf - 2f. 

As a further check that the fluxon lattice really moves 
as a unit, we have used eq. to compare the time- 
averaged voltage (V) at I/I c = 0.04, on the upward cur- 
rent sweep just above the first voltage jump, to the di- 
mensionless period At of the a. c. component. The pre- 
diction for (V)/(RI C ) from eq. 10 (namely (V)/(RI C ) = 
0.0381) is in excellent agreement with the value calcu- 
lated by a direct time-average of (V)/(RI C ) (0.0372). 
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FIG. 8: (a) Normalized flux vs. plaquette number for the 
parameters of Fig. [7] (b) Ground state vortex charge con- 
figuration for the same parameters as (a), Only 40 of the 120 
plaquettes are shown. 




FIG. 9: IV curve for / = 41/120, Qj = 10.0 and Aj = 
1.0, and N = 120. There are three jumps on the increasing- 
current branch of the IV characteristics; the first occurs near 
I/I c ~ 0.01, and the others can be easily seen in the Figure. 
The inset shows the IV curve for < I/I c < 0.05. The 
interpretation of these regimes is given in the text. 



E. f = 1/3 +1/N 

Next, we consider the behavior of the same system 
with one additional fluxon. Fig. shows the IV curve, 
and Fig. ITUT a) displays the normalized flux vs. plaque- 
tte number at I/I c = 0, for / = 1/3 + 1/120. In this 
case, the additional fluxon fractionalizes into three exci- 
tations. These are visible in Fig. HOf a) as regions which 
are distorted, in comparison to the flux plot at / = 1/3. 



As in Fig. Eta), these regions show how the / = 1/3 
flux pattern changes to accommodate the additional flux 
quantum. To understand why this extra fluxon fraction- 
alizes into three parts, we show in Figs. llOf b) and (c) two 
possible patterns of vortex number when the additional 
fluxon is added. The / = 1/3 ground state consists of a 
+1 vortex every three plaquettes. If the additional fluxon 
is un-fractionalized, we expect that it will enter the lad- 
der as in Fig.^|(b). This pattern must have two neigh- 
boring plaquettes with a +1 vortex number, and hence is 
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FIG. 10: (a) Normalized flux vs. plaquette number for the 
same parameters as in Fig. Q except that / = 41/120. (b) A 
possible high-energy vortex charge configuration for this field, 
(c) Schematic of the actual ground-state vortex charge con- 
figuration, as obtained after iterating the equations of motion 
up to a time T max = 5000 and I/I c = 0. In both (b) and 
(c), we show only 40 of the 120 plaquettes. In the actual lad- 
der at / = 41/120, the three 1/3-fluxons are separated by 40 
plaquettes, as suggested in Fig. 10(a). 



a relatively high-energy state. It is much more plausible 
that the lowest-energy pattern would never have +1 vor- 
tices in adjacent plaquettes. Such a pattern is shown in 
Fig. [23 (c). The 1/3 fluxon shows up as a region where 
two +1 vortices are separated by only one empty pla- 
quette, rather than two as in the unperturbed / = 1/3 
pattern. To conserve the number of +1 plaquettes, there 
must be three such regimes; hence, the added fluxon is ex- 
pected to fractionalize into three 1/3-fluxons, as we find 
numerically. These 1/3-fluxons are equidistant in order 
to minimize the repulsive energy between these added 
excitations. 

When the extra fractionalized fluxon moves through 
the lattice, it generates voltage. Once again, we can use 
the analog of eq. (|11|) to connect the time-averaged volt- 
age to the oscillation period of V(r). This analog is 



(V) _ 3 ■ (2tQ 
RI C NQjAt' 



(12) 



Here At is now the period of the a. c. voltage in the 
regime where three 1 /3-fluxons are moving, measured in 
units of 1/ojp. When we use this relation to calculate (V) 
from our calculated At, at I/I c = 0.02 (in the regime of 



three moving 1/3-fluxons), we find that the value of (V) 
obtained from eq. (|12[1. using the calculated value of Ar 
is in excellent agreement with the value of (V) calculated 
directly. 



F. f = 1/q and f = 1/q + 1/N, with q > 3 

We have also carried out similar calculations for 1/4, 
1/5 and 1/6. The ground state (at I/I c = 0) consists, 
as expected, of a periodic array of +1 vortices separated 
by q — 1 zero vortices (q = 4, 5 or 6). When we carry 
out the corresponding calculations for / = 1/4 + 1/iV, 
/ = 1/5 + 1/N and 1/6 + 1/N, we again find that the 
added vortex fractionalizes into q equally spaced pieces, 
as suggested by our results for q = 2, 3. 



G. Effects of Varying Ladder Length, Inductive 
Coupling, and Damping Coefficient. Effects of 
Removing One Fluxon. 

We have also studied / = 1/q + 1/N (q > 1) for a 
ladder of double the length, i. e., 240 plaquettes. The 
flux patterns are very similar to those of the shorter lad- 
ders, except that, of course, the fractionalized fluxons are 
twice as far apart as for N = 120; the IV curves are also 
qualitatively similar. 

As another important test case, we have computed 
both the ground state configurations and IV character- 
istics for / — > / — 1/N, in order to compare the effects 
of removing rather than adding one fluxon. (We do not 
display these results graphically.) In general, the results 
are similar to those found for adding one fluxon. The 
results for / = 1/2 — 1/N are nearly identical to those 
for / = 1/2 + 1/N: the missing fluxon is fractionalized 
into two half-fluxons. For / = 1/q — 1/N, with q > 1 
an integer, the IV and flux characteristics are similar to 
those for / = 1/q + 1/N; however, the flux curves are 
not identical, except for q = 2. It is not surprising that 
one obtains different results for adding and subtracting 
a fluxon for q > 2: only at q = 2 would one expect, on 
symmetry grounds, that adding and subtracting a vortex 
would give the same response. 

We have also tried varying the strength of the inductive 
coupling constant Xj, considering Aj = 0.01, 0.1 and 10.0 
as well as 1 for / = 1/2. For Aj = 0.01, the system 
behaves like a collection of independent junctions: (V) — 
for all I/I c < 0.9, and the flux pattern is very similar to 



that of Fig. A j = 0.1 also gives results very similar to 
those of Fig. 21 But for A j = 10, the fluxons are depinned 
even at I / I c = 0, and the IV curve is very similar to that 
seen in Ref. |26j for a long junction in the continuum 
limit. We have also carried out calculations (not shown 
here) for values of Qj = 1.0, 5.0 and 50.0 for / = 1/2, 
N = 120 and Aj = 1.0. The position of the first jump 
in the IV curve remains unchanged even for the critically 
damped case (Qj = 1.0). We found that, for the larger 
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FIG. 11: Normalized flux vs. plaquette number, plotted for 
N = 121 and / = 61/121. 
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FIG. 12: IV characteristics for N = 121 and / = 61/121, for 
the parameters of the preceding figure. Inset: enlargement of 
IV characteristics at small currents. 



Qj's, the fluxons are expelled from the ladder at lower 
currents than for smaller Qj's. Finally, for Qj = 1.0, the 
IV curve of the ladder closely resembles that of a single 
small, critically damped junction. 



H. f = 1/2 + 1/(2N): One Added Half-Fluxon 

As a final calculation, we have considered a ladder 
with an odd number of plaquettes, namely N = 121, 
and / = 61/121. Intuitively, one might expect that this 
ladder would correspond to a single added 1/2-fluxon su- 
perimposed on an / = 1/2 ground state. The latter 
corresponds to / = 1/2 = 60.5/121, while the extra 1/2- 
fluxon is needed in order to bring the total flux up to 
/ = 61/121. Indeed, the flux plot, shown in Fig. ^] is 
consistent with this picture: an alternating flux pattern 
of the / = 1/2 ground state on which is superimposed 
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> 
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4020 



FIG. 13: Time-dependent voltage V(t)/(RI c ) plotted for di- 
mensionless time r between 4000 and 4020, for f = 61/121 and 
current I/I c = 0.03. Qj = 10 and A 2 7 = 1. The dimensionless 
period of oscillation At satisfies eq. 1131 

)• 



the flux pattern of a single 1/2-fluxon, i. e., one of the 
two excitations distinguishable in Fig. 6. 

The IV characteristics of this ladder (shown in Fig. I12|) . 
are consistent with this picture. There are three regimes 
of non-zero voltage (seen most clearly on the decreas- 
ing current branch), namely (i) motion of the 1/2-fluxon 
through the background of the pinned / = 1/2 lattice; 
(ii) motion of the depinned f — 1/2 lattice; and (hi) 
purely resistive behavior, with all small junctions in the 
voltage state. The motion of the half-fluxon [in regime 
(i)] is particularly intriguing. At a given current, the volt- 
age generated by this 1/2-fluxon is ~ 1/2 that produced 
by the two moving 1/2-fLuxons at the corresponding cur- 
rent in Fig. 6. This result is consistent with the picture 
that the two 1/2-fluxons in Fig. 6 move at the same av- 
erage velocity at any given applied current. In addition, 
as mentioned earlier, at all currents in this regime, the 
time-averaged voltage (V) generated by the 1/2-fluxon is 
related to the ac frequency v of that voltage by 



(V) _ 2ir 
W c ~ NQjAt' 



(13) 



where At = oj p /u. As an illustration, we show in Fig. 1131 
V{t)/{RI c ) for / = 61/121 and a ladder of length 121 
plaquettes, at a current I / I c = 0.03, in the regime where 
the voltage is produced by a single 1/2-fluxon. 



IV. DISCUSSION 

We have shown that a single fluxon introduced into a 
Josephson ladder will fractionalize if the ladder is initially 
in an / = 1/q ground state with q = 2, 3, 4, 5, or 6. This 
fractionalization leads to a characteristic behavior of the 
(suitably defined) flux in the ground state at I/I c = 0, 
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and even for I/I c ^ 0, provided that I is smaller than 
the critical current for onset of voltage. Once a measur- 
able voltage is obtained across the ladder, for these fields, 
there is a characteristic series of regimes with increasing 
current: first, the fractionalizcd fluxon moves, generating 
a voltage; then, the underlying fluxon lattice is depinncd; 
and finally, the ladder as a whole switches to the voltage 
state and the fluxons are expelled from the ladder. 

The type of fractional excitation found here are ba- 
sically walls between different domains of the underly- 
ing 1/q ground state. For an applied field of 1/q, there 
are q discrete flux lattice ground states, each shifted by 
one plaquette from the next ground state. If two dif- 
ferent ground state domains, shifted by one plaquette, 
are placed next to one another, a 1/q fractional fluxon 
is formed. Such walls, with fractional vortex charge, 
are the analogs of fractional excitations well known in 
many other physics contexts, especially in one dimen- 
sion. For example, solitons in polyacetylene are frac- 
tionally charged excitations which can form between two 
domains of alternating singly and doubly bonded car- 
bons in this one-dimensional chain compound 27]. The 
fractional fluxons in the Josephson ladder are unusual 
because it appears that a fluxon of any rational fraction 
can form in a suitable applied magnetic field. 

The fractional vortices described in this paper gener- 
ally produce a non-hysteretic IV characteristic, even if 
the individual junctions which make up the ladder are 
hysteretic. This behavior is similar to that previously ob- 
served for integer vortices in two-dimensional arrays psj. 
which also move in a non-hysteretic fashion even when 
the junctions in the array are hysteretic. In the case of 
the 2D arrays, the explanation for this behavior is that 
the effective quality factor describing the voltage gen- 
erated by the moving vortex is Q e ff = {2eR 2 IdC /K) 1 / 2 , 
where C is the junction capacitance, R the junction shunt 
resistance, and Id is the current required to depin a vor- 
tex from its equilibrium position within a plaquette of 
the 2D lattice. Since Id ~ 0.1/ c in the 2D case, where 
I c is the critical current of one junction, Q e ff is much 
smaller than the junction quality factor (2eR 2 I c C/H) 1 ^ 2 , 
and hence the motion of the vortex is overdamped even 
when the individual junctions are underdamped and hys- 
teretic. 

For the ladder arrays, the depinning current Id for a 
fractional fluxon is less than 0.017c for all cases we have 
considered in this paper. Hence, even though Qj = 10, 
the analogous effective quality factor for the fractional 
fluxon will be less than 0.1. The resulting IV character- 
istics should thus be non-hysteretic and similar to those 
of an overdamped Josephson junction, as we have ob- 
served numerically (see Fig. 13). However, we have not 
succeeded in deriving a driven-pendulum-like equation of 
motion for the fractional vortex, similar to that given in 
Refs. 25] and 28]. 

One striking feature of these fractionalized fluxons is 
that they do not appear to form in two dimensions. When 
a fluxon is added to an f — 1/2 ground state of a square 



lattice in two dimensions, the fluxon is depinned at a 
current much lower than that of the / = 1/2 lattice, just 
as in one dimension. However, numerical studies suggest 
that the fluxon is not fractionalized, but remains a com- 
pact object 28]. Fractionalization is apparently easier in 
ID because the domain walls are zero-dimensional, and 
hence require much lower energy to form. 

It is of interest to compare this prediction to what we 
would expect from a single un-fractionalized fluxon in an 
otherwise empty ladder. Here, the expected mechanism 
of fluxon motion is that the fluxon moves by one pla- 
quette. An argument similar to that just given yields 
(V) = [h/{2eN)v. 

We can also consider the motion of an un-fractionalized 
vortex through the background of a pinned / = 1/q 
fluxon lattice. Since the motion of the single vortex 
would have a period of q lattice constants, and since 
there would be a phase slip of 2ir each time the vor- 
tex circulates once around the ladder, the equation re- 
lating the time-averaged voltage to the period would be 
(V)/(RI C ) = 27rq/(QjAr), where At is the period of the 
ac voltage, in units of 1/lu p . This relation is the same as 
would hold for a vortex fractionalized into q parts and 
moving through the same background. 

Thus, the clearest experimental way to see a signature 
of the fractionalization would be to compare the period 
of the ac voltages with the time-averaged voltage in a lad- 
der with isolated fractional vortices. Such ladders were 
discussed earlier - an example is the case of N = 121 and 
/ = 61/121, which contains a single 1/2-fluxon in a back- 
ground of an / = 1/2 fluxon lattice. Besides searching 
for a signal of the fractionalized vortices in the IV char- 
acteristics, perhaps the flux profile of the fractionalized 
vortices could also be directly imaged. 

As mentioned earlier, the fractional vortices discussed 
here were previously suggested ^| for a ladder of alter- 
nating and tt junctions. Since fractional vortices can 
thus appear in both types of ladders, one might ask if 
any information about the nature of the junctions (i. e., 
whether they are or tt junctions) could be gleaned from 
a simple inspection of the IV characteristics. If all the 
rungs of the ladder are tt junctions, they would have the 
same IV characteristics as a ladder made entirely of 
junctions, since the presence of tt junctions would pro- 
duce no additional frustration in this case. Thus, the IV 
characteristics would not distinguish between all and 
all 7r junctions. If the ladder contained alternating and 
tt junctions, the IV characteristics, in the presence of a 
magnetic field, would differ from that for all (or all tt) 
junctions by being shifted in magnetic field by 1/2 a flux 
quantum per plaquette - that is, the IV characteristics of 
the alternating ladder at field / would be the same as that 
of the or tt ladder at / + 1/2. For a more complicated 
distribution of and tt junctions, the IV characteristics 
of the ladder would be more complicated, and the pres- 
ence of fractional vortices might be difficult to detect. 
Thus, the existence of fractional vortices, of itself, might 
not be sufficient to determine if the ladder were partly 



10 



composed of ir junctions. 

Finally, we note that the present results are all ob- 
tained from a classical set of governing equations. If 
the individual junctions are sufficiently small, we expect 
that the fractional fluxons would behave like quantum- 
mechanical particles. Specifically, the phase config- 
uration of the ladder would be described by an N- 
dimensional wave function ^((fii, ....</>jv), which would be 
an eigenfunction of the ladder Hamiltonian (expressed as 
an operator). It would be of interest to calculate the spec- 
trum of such eigenstates, and to study the behavior of a 
lattice, in the small junction regime, containing two or 
more such quantum-mechanical fractional fluxons. One 
important question would be the statistics (bose, fermi, 
or fractional) obeyed by such particles. Such questions, 



though fascinating, are, however, beyond the scope of the 
present paper. In addition, there should be a broad range 
of experimental parameters where the classical equations 
are applicable. Thus, the present predictions should be 
readily testable experimentally. 
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